Generalized synchronization of chaos in non-invertible maps 
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The properties of functional relation between a non-invertible chaotic drive and a response map 
in the regime of generalized synchronization of chaos are studied. It is shown that despite a very 
fuzzy image of the relation between the current states of the maps, the functional relation becomes 
apparent when a sufficient interval of driving trajectory is taken into account. This paper develops 
a theoretical framework of such functional relation and illustrates the main theoretical conclusions 
using numerical simulations. 
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I. INTRODUCTION 

Since the first studies by Van der Pol the dynamical 
theory of forced synchronization relates the synchroniza- 
tion phenomena to the onset of stable response behavior 
of driven oscillator Jj]-[|. As a result of this stability, 
the beats in the response oscillations disappear and a 
stable periodic motion occurs in the phase space of the 
driven oscillator. After that the oscillator becomes en- 
slaved by the periodic forcing. Two main bifurcation 
scenarios that lead to the formation of the stable re- 
sponse behavior are the birth of stable limit cycle on 
a torus and Andronov-Hopf bifurcation |0,@,0|. 

Based on these theoretical frameworks recent studies of 
synchronization in the chaotic oscillators have lead to the 
development of various notions of chaos synchronization. 
These notions include: identical synchronization [fTO fl2| |, 
generalized synchronization 13-1(|, and phase synchro- 
nization 

An extension of the stability approach towards the case 
of directionally coupled non-identical chaotic oscillators 
is captured by the notion of generalized synchronization. 
The term, generalized synchronization, was introduced 
in pj and used to describe the onset of synchronization 
in directionally coupled chaotic systems as the formation 
of a continuous mapping that transforms a trajectory on 
the attractor of the drive system into a trajectory of the 
response system. For the systems with invertible dynam- 
ics this is equivalent to the formation of a continuous 
mapping which links the current states of the systems 
when they are settled down on the synchronous attrac- 
tor. 

In the case of invertible dynamics of the driving system 
the relation between attracting properties of the response 
behavior and the some of the properties of the synchro- 
nization mapping has been reviled and proved |l9|-|2~4j| . In 
the case of non-invertible dynamics the relation between 



the response stability and existence of the synchroniza- 
tion mapping was indirectly observed with the auxiliary 
method J25j,^(|. However, the detailed and rigorous study 
of this relation has not been done. 

This paper presents the results of theoretical and nu- 
merical study of the onset of functional relation between 
the chaotic trajectories of driving non-invertible system 
and the states of the response system. Some results about 
this type of synchronization can also be found in p7| , ^8| . 
To be specific we consider the systems in the form of the 
following maps 



/(•En)} 
Vn+1 = 9c(Xn,y n ), 



(1) 

(2) 



where equation (Q) describes non-invertible driving sys- 
tem and equation (|J) the response system. We assume 
for the sake of dcfinitcness, that in the system (Q), (Q) 
x e X C 5R m and y e Y C 5R £ , and that / and g c are 
continuous. The subscript c stands for a coupling param- 
eter. The global dynamics generated by the system (Q), 
(0) can be written in the form F c (x n ,y n ) = (x n+1 ,y n+1 ). 



II. THEORETICAL RESULTS 

The non-invetibility of / implies that for any xq 6 
X, there are possibly infinitely many different backward 
orbits of xq (xq, X-i, ...). Let us denote by P(xq) 

the union of all these possible backward orbits. Denote 
by p one of these orbits. 

Assume that there exists a partition {Xi, Xi} of X, 
such that f\xj is one-to-one to the imageQ We define 
the backward symbolic sequence a :— (cto, cti, a n , ...) 



'The expression f\xj stands for the restriction of function / 
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associated to a backward sequence (xo, X—i, X- n , ...) 
with 014 = j if x-i S Xj. 

Since we study synchronization in dissipative systems 
we also assume that there exists a ball of dissipation 
B C W n+e , i.e. F C (B) C Int(B) for any c e 5, where 
S" is a region in the coupling parameter space in which 
system (^) has stable response behavior. Without loss 
of generality we assume that B = B x x B y , i.e. B is a 
rectangle, where B x (resp. -B y ) is a ball in the x-space 
(resp. y-space)). Denote by A c the maximal attractor in 
B, i.e. A c = r\^ =0 F c( B )- Assuming that system (f| has 
stable response behavior, we have 



lim \y„ 



Vn 



0. 



(3) 



where (x n ,y n ) = F?(x ,y ) ((x n ,y n ) = F?(x ,y )) and 
{xo,Uo), (xo,Vo) are arbitrary points in B. 

Let A c , x ■— n x A c be the image of A c under the natural 
projection II X to X. The set A c , y is the image of A c by 
IL,, the natural projection to Y. 

Theorem 1. Under assumption of stable response 
given by Eq.(^ the attractor A c is the union of graphs 
of infinitely many functions. Each function h a is deter- 
mined by a symbolic backward orbit a. Moreover each h a 
is continuous. 

Scheme of the proof. The main point of the proof is 
to note that Xq and a determine all the backward orbit 
X— i, ...,x_ n , ...). From h ere the proof is similar 
to the one for the invertible case [^| . Given p, one can 
define h a (xo) as the following limit: 

h a (x ) := lim Tl y F?(x^ n ,y) 

n — >oc 

independently of y G B y . This limit exists because of 
assumption 

Continuity is proved in the same way that in the in- 
vertible case taking into account that a is the same for 
x and x, two close points in X (see 23 for details). □ 



Assume that the driving system is generated by a map 
/ of the interval /. Then we can say more about regular- 
ity of branches provided that \y n — y n \ goes to zero ex- 
ponentially fast. Indeed the following proposition holds. 

Proposition 1. Assume that 



(4) 



where A > and < A < 1 are constants depend- 
ing on c. Assume also that the one- dimensional driv- 
ing system f is "hyperbolic", i.e. there exists uq such 
that (f n )'(x) > -p- > 1 for all n > no and for any x 



for which the derivative exists. Then each branch h a is 
Lipschitz continuous. 
Scheme of the proof. 

The proof is basically the same of the one of theorem 
7 in [ p3[ . The point is to replace the space of Lipschitz 
functions £/l,m by Hl^m,o with the restriction that the 
domain for h a € Hl^m,u consist of values of x compatible 
with a. 

Remark that by assumption, one has 7_ < 1. The 
condition (16) in theorem 7 in [E3| can be replaced by 



< AX n < 



1 



7- 

which is always true if n is large enough. □ 



(5) 



Taking into account details of the proof of this propo- 
sition (which are omitted here and can be found else- 
where [g3[) one might expect that if / is not hyperbolic 
then continuous synchronization functions h a could be 
not the Lipschitz-continuous functions. In this case the 
picture of individual branches h a might contain wrinkles 
and cups, and as the result, might appear fuzzy in nu- 
merical simulation due to finite resolution. 

From now on we denote by h the union of all graphs 
h a . Each branch h a is only defined for the values of 
x compatible with a. That is, if there is no backward 
orbit for the value x with the symbolic itinerary a, then 
h a is not defined at x. Thus, in general, the number of 
branches of h is not necessarily the same for all x. 

Corollary to Theorem 1 If {Xi} is a Markov par- 
tition then for any admissible a, the domain of h a con- 
tains an element of {Xi}, i.e. if a — (cxq, a n , ...) then 
D{h a )DX ao . 

One may ask when two different branches of h are close. 
The next theorem shows that two branches that have sim- 
ilar recent history, are close to each other. Branches with 
different recent symbols could be either close or far from 



each other. As numerical simulations discussed in Sec. Ill 
we show that branches with different recent symbols may 
even intersect each other. 

Theorem 2. Let {a 1 } be a series of infinitely long 
symbolic sequences a 1 £ {1, Z} N such that lim^oo a 1 — 
a in the standard product topolog 1 ^ Jjjj/ . Then, 

lim \h y - h a \ = , 

where \h fi \ := sup xeD{p) \h^(x)\ 
Proof. 



to the partition element Xj . This type of partition is standard 
in the Markov maps of the interval, see for example |J2Sfj . 



^This means in particular that the greater is i the greater TV 
first symbols in a and a 1 coincide. 
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From (g), there exists an N such that for n > N and 
for any y and y 



\y n -y n \<e 



(6) 



This is, of course, true also for yo, yo £ Ac, y - 

Now, choose i big enough in such a way that the 
first iV symbols of a and a? coincide. This means that 
(x, h a {x)) and (x,h a (x)) have the same first iV x-pre- 
images (x_i,...,x_at). Then one can apply (^) for these 
points obtaining 

\h &i (x)-h a (x)\ < e . 

This is true for any x compatible with a and a 1 . The 
theorem is proven. □ 

If one assumes monotonicity in assumption (|^), that is, 
\y n — y n \ < \Y\c n , then a stronger result holds: Let a and 
a have the same first n symbols, then \h a — h a \ < \Y\c n . 

Let us remark that there are three logical possibilities 
for the structure of the set h(x): it could be finite, it 
could be countable, and it could contain a Cantor set. 

The first possibility trivially occurs in the case of 
identical synchronization, that can exist even for non- 
invertible / (see for instance p3|). 

The second possibility can be justified by construc- 
tion of a special pair /, g c . We believe that it could be 
done. It is definitely true for the case when / has zero 
topological entropy, i.e. the number of admissible words 
{ao, ...a„_i} grows subexponentially as n goes to infin- 
ity. For example, if the topological Markov chain corre- 
sponding to a Markov partition {Xj} has zero topological 
entropy then the number of admissible words grows not 
faster than polynomially. If topological entropy is posi- 
tive then the construction of such an example is rather 
difficult, we are going to study this problem elsewhere. 

The third possibility seems to be generic in a space of 
pairs of functions /, g c provided that / has positive topo- 
logical entropy. It seems natural to believe that differ- 
ent driving signals correspond to different outputs in the 
regime of synchronization. Numerical simulations con- 
sidered in Sec. Ill confirm this conjecture (see also Fig. 3 
in |27| and Fig. 2d in |p8[). However, we do not have a 
rigorous proof of this statement right now. 

The statements made above give a clear picture of 
the underlying structure of generalized synchronization 
of chaos in the case when the dynamics of drive system 
is non-invertible. 



of logistic map driven by tent map. In this case function 
/(x„) in the driving system (pi) is of the form 



x/b if x < b, 

(x-l)/(b-l) ifx>6, 



where 6 is a control parameter, < b < 1. We will con- 
sider the case when the dynamics of the response system 
is given by the following map 



Vn+i = (1 - e)ay n (l - y n ) + e/(x n ), 



(7) 



where a is the control parameter of the map and e £ 
[0, 1] is a coupling parameter. Note that upper bound 
of the contraction rate in the y-direction denoted by c 
in previous section is, in this case, c = (1 — e)a. In the 
numerical simulation considered in this section the values 
of control parameters are fixed b = 0.677 and a = 3.7. 

Since the dynamics of the maps are different, the gen- 
eralized synchronization is the only possible regime of 
synchronization, except to the trivial case when e = 1. 
The onset of generalized synchronization is detected with 
auxiliary system approach. In this analysis we study the 
stability of the chaotic response attractor in the mani- 



fold y n 



where variable z n is described by an exact 



replica of the system 
tern [E5|. 



0.0 



which is called auxiliary sys- 
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III. SYNCHRONIZATION FUNCTION IN 
NON-INVERTIBLE MAPS: EXAMPLE 

To illustrate the properties of generalized synchroniza- 
tion in non-invertible maps consider the synchronization 



FIG. 1. The dependance of conditional Lyapunov exponent 
A on the value of coupling parameter e - (a), and the values 
of maximal and rms deviations of d n computed as a function 
of e - (b). 
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It can be shown with the analysis based on contraction 
mapping ||(J that stable identical oscillations in response 
and auxiliary systems are guaranteed when the values of 
coupling parameter are within the interval 1 — 1 / ai < e < 
1 + l/a r , where a t = a(l - 2y™ m ), a r = a(l - 2y™ ax ), 
and y™ m and y™ ax are the leftmost and rightmost points 
of the attractor in the response system. This indicates 
that if for the selected parameter values of the maps the 
coupling is stronger than e = 0.7, then synchronization 
is monotonically stable. 

More precise evaluation of the synchronization thresh- 
old can be done with the analysis of conditional Lya- 
punov exponents (see, Fig |l|a) and with the analysis 
of deviation of response-auxiliary system from manifold 
Vn = z n- The dependance of maximal, d max , and rms, 
drms, values of deviation d n = z n — y n on the values of 
coupling parameter are presented in Fig. [j]b. 

Based upon the plots in Fig. |l| one would expect that 
the regime of generalized synchronization takes place for 
the coupling parameter values e > 0.3 This regime as- 
sumes existence of a continuous functional relation be- 
tween the trajectories x n and y n . However, when one 
plots the states of the one-dimensional phase space y n 
on the synchronized attractor versus the corresponding 
states of one-dimensional phase space x n , the image of 
such relation looks fuzzy, see Fig. 0. It is clear that point 
to point mapping x n — » y n cannot be interpreted as a 
continuous function. 




FIG. 2. Synchronized chaotic attractor computed for 
e — 0.6 plotted in the phase plane (x n ,y n ). 

In order to unveil the synchronization function 
in the numerical analysis we consider the mapping 
(x n , [a\, ...,a m ]) — > y n , where are symbolic represen- 
tation of prehistory of the trajectory x n . These symbols 
are generated by the tent map (|l|) which is partitioned 
into two regions: = L if i„_t < b and Qfe = R 



if 



> b. The points of the synchronous attractor 



whose preceding symbolic sequence of the driving tra- 
jectory has a specific sequence (a mask) were selected 



and analyzed separately. This analysis shows that as 
the length of the sequence increases the cloud of points 
shrinks into a curve. An example of such convergence 
is illustrated in Fig. |], where the mask of 8 symbols 
[L, L, R, L, L, R, R, R] is studied. Similar behavior is ob- 
served for the other sequences of the same length. In the 
sequences of this length the shape of the curve varies as 
the symbolic sequence changes. All together these curves 
form the fuzzy shape of the synchronized attractor shown 
in Fig.| 
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FIG. 3. Points of the synchronized chaotic attractor shown 
in Fig. ^ generated by the trajectories which symbolic se- 
quence of length m — 8, preceding the final point x n , fits to 
the mask [L, L, R, L, L, R,R,R]. e = 0.6 

To evaluate the convergence of the image of synchro- 
nized attractor to a continuous function h a (x) we ana- 
lyzed the sets of attractor points conditioned by all possi- 
ble symbolic masks a of various length m. For each mask 
of preceding symbols S l m = [a\, .., a m ] we computed the 
best polynomial fitting function (pgi (x) of order 6 us- 
ing Singular Value Decomposition (SVD) algorithm, and 
studied the dependance of mean squared error (MSE), 
averaged over all masks of length m, versus m. This de- 
pendence computed for the values of coupling parameter 
e = 0.6 and e = 0.4 are shown in Fig. ^. 

One can see from the Fig. || that MSE decreases ex- 
ponentially fast when m increases. Approximating this 
dependence with exponent 



MSE(m) 



(8) 



one can find the rate of convergence, which is in the case 
e = 0.6 equals A w -1.05. 

Figure shows how the convergence rate A depends 
on the value of coupling parameter e. Comparing this 
plot with the plot of conditional Lyapunov exponent ver- 
sus e one can clearly see the similarity in these plots. 
This is indicative of the fact that in the generalized syn- 
chronization regime in our case the convergence rate A is 



4 



related to the contraction rate of response system given 
by conditional Lyapunov exponent A. This relation was 
defined for the case of monotonic stability of the response 
behavior, see Section O. 
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FIG. 4. The dependence of MSE of best polynomial fitting 
function for the attractor points y n ,x n on the length m of the 
preceding masks 



explains the formation of fractal structure in the syn- 
chronization mapping which consists of infinite number 
of the branches. 
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FIG. 6. Branches of the synchronized chaotic attractor 
shown in Fig. ^| computed for symbolic sequences of length 
m = 8, [on , «2, R, L, L, R, R, R\. The values of three most re- 
cent symbols ai,a2,R are shown next to the corresponding 
branch. 
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FIG. 5. The dependence of the convergence rate A on the 
value of coupling parameter e. 

One can see from Fig.^ that synchronization mapping 
for the current states of the drive and response system 
has a complex structure of branches. The appearance of 
Cantor set in the synchronized attractor caused by non- 
invertible driving system was reported before in |2|] . The 
formation of the fractal structure can be explained by the 
analysis of deviations of the branch caused by the change 
of symbol a n appeared in the symbolic sequence n itera- 
tions before. For the most of the driving trajectories, the 
stability of the response behavior acts in a such way that 
the deviation will be reduced with the increasing value 
of n. As the result the deviations of different scales for 
different values of n are responsible the formation of the 
complete structure of the synchronized attractor. This 
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FIG. 7. Branches of the synchronized attractor for the se- 
quences [L, ct2, «3, L, L, R, R, R\. The values of three most 
recent symbols L, ot2, 03 are shown next to the corresponding 
branch. 

This mechanics behind the formation of Cantor set is 
illustrated in Figs. || and |?]. Figure || shows the set of 
points of the synchronized attractor corresponding to the 
trajectories whose 8 most recent symbols are fixed. One 
can see that branches with the same most recent symbols 
remain relatively close to each other, while the branches 
with different most recent symbols can get apart signif- 
icantly. The similar situation takes place for the alter- 
nation of the more remote symbols, see Fig. |^. How- 
ever the scale of the deviation in this case is less then in 
Fig. ^. Further increase of the length of alternating sym- 
bols by one in considered mask of symbols the number 
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of branches doubles. The new branches appease close to 
the branches with the same most recent symbols. Con- 
tinuation of this process with infinitely long mask n — > oo 
leads to the formation of a fractal, Cantor type set. 



IV. CONCLUSIONS 

The results of the theoretical analysis of synchroniza- 
tion presented in Sec. O and the numerical analysis of 



particular example considered in Sec. Ill allow one to 



draw a number of important conclusions on typical prop- 
erties of synchronization mapping that characterize gen- 
eralized synchronization of chaos in the case of non- 
invertible driving system. 

The synchronization mapping, in this case, can be in- 
terpreted as a continuous function only when all back- 
ward iterations X-„ of the driving trajectory xo are in- 
cluded in the vector of the function arguments. We have 
shown that Markov partition of the driving map can be 
used to describe the backward iterations in a compact 
symbolic way. The set of points of synchronization map- 
ping plotted for the trajectories, which prehistory is con- 
ditioned by a selected sequence of symbols, asymptoti- 
cally approach the graph of a continuous function as the 
length of the sequence increases. 

The set of these graphs given by all possible symbolic 
sequences of infinite length form the complex fuzzy object 
which is typically observed in the joint phase space of 
the one-dimensional chaotic maps synchronized in the 
generalized sense. 

An interesting feature of the considered example with 
1-d maps is that branches h a remain smooth functions 
even when the rate of conditional stability is rather weak. 
Indeed, according to the papers |]l9|-^| , one would expect 
that h a would become non-differentiable, Holder contin- 
uous function as the contraction rate in the response sys- 
tem becomes lower than some critical value. However, 
this critical value is given by a contraction rate towards 
the chaotic attractor in the driving system. In our case 
the driving trajectories of the 1-d map with the specified 
symbolic sequence a do not have contracting direction. 
They have only unstable direction. As the result the con- 
traction rate in the response system for any given branch 
h is always stronger than in the driving system. 

This feature is typical for the driving systems in the 
form of one-dimensional hyperbolic map. In the case 
of non-hyperbolic 1-d driving map function h a may be- 
come non-differcntiable and contain wrinkles and cusps. 
Another interesting question arises when one considers 
a system without exponential stability or Markov parti- 
tion? It could be expected that branches become non- 
differentiable. For the moment the question is open. 
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